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(54) Abstract Title: Calculating an estimate of bit reliability In a OFDM receWer by multiplication of the channel 
state modulus 



(57) A scheme for simplifying the computational 

complexity of calculating log likelihood ratios (LLR) 
for soft output demapping is provided. The scheme 
can be implemented inside a receiver in an OFDM 
system using multi-level modulation whereby 
calculation of LLR ratios is accomplished using only 
multiplication of the channel state (CSI) modulus by 

(which represents the half distance between 
partition boundaries in a QAM constellation). 
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At least one drawing originally filed was informal and the print reproduced here is taken from a later filed formal copy. 
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A METHOD OF CALCULATING AN ESTIMATE OF RFXl ABILITY 
OF BIT INFORMATION IN A OFDM RECEIVER AND A RECEIVER 
OPERATING IN ACCORDANCE WITH THE METHOD 

Backeround to the Invention 

The present invention relates to a method of estimating a reliability measurement for a bit 
to be decoded in a multi-level transmission scheme based on bit-interleaved coded 
orthogonal frequency division multiplexing. The method significantly reduces the 
computational task of the receiver. The present invention also relates to a receiver and 
demapper operating in accordance with the method. 

A significant amount of effort has been expended in the development of high data rate 
digital communication. Such communication technologies have brought about wireless 
local area networks, such as those defined in IEEE802.1 la and HIPERLAN/2. 

Reducing the computational complexity in the receivers of such a network would be 
beneficial for implementation of such systems because of the limited resources available in 
many low cost commercial devices. 

Figure 1 schematically illustrates a communication system wherein the present invention 
can be applied. The transmitter part is in accordance with the HIPERLAN/2 and 802.1 la 
standards. In both of these standards difierent QAM modulation formats are used, such as 
BPSK, QPSK, 16 QAM and 64 QAM. 

In general a first user (who may be human or machine) wishes to transmit a data stream to a 
second user. The bits in the data stream are firstly convolutionally encoded by an encoder 
2 to implement forward error correction. Schemes for convolutional encoding are well 
known to the person skilled in the art and hence only a brief description is required here. In 
general convolutional encoding is perfonned by shifting the data through a chain of delay 
elements usually implemented as a shift register. In the case of the HlPERLAN/2 and 
802.11a standards the convolutional code has a rate Yi and a constraint length of seven. 
The chain is, as shown in Figure 2 tapped at various points along it and the tapped signals 



2 

are summed together at summers 4 and 6 which are implemented as modulo 2 adders, 
which are often physically realised as a cascade of exclusive - OR gates, 

« 

The arrangement shown in Figure 2 has six delay elements T and data can also be tapped 
off at the input to the first element. The encoder shown has, in effect, seven snapshots of 
the input string and is therefore described as having a constraint length of 7. 

The encoder produces two outputs Yi and Y2 in accordance with the series of connections 
and taps that are formed within it. Working from the input to the output, the first output is 
represented by the first, the third, fourth, sixth and seventh taps, or 101 1011 . The second 
output Yi is represented by 1 1 1 1001. Thus the code generator for this encoder is [133, 
171] in octal format. It should be noted that since the convolutional code has rate '/i the 
data entering the encoder is a scalar bit stream, whilst the data leaving the encoder is a 
stream of binary bit pairs. The data ou^ut of the encoder 2 is sent to a bit interleaver 10 so 
as to spread the bite so that potential errors caused by channel degradation are seen as 
independent at the receiver. The interleaved bits are then converted into quadrature - 
ampUtude modulated (QAM) signals at a QAM m^r 12. 

In other words sequences of 2m bits are assigned to one of the M=2^ points in the QAM 
constellation. The mapping is performed as a Gray coding such that sequences of bite 
associated with adjacent ^robols only differ by one bit. 

An exemplary 16-QAM constellation is shown in Figure 3. 

The complex (real and imaginary parte) signals are then fed to an orthogonal frequency 
division multiplexing (OFDM) modulator 14. OFDM is a multi-carrier modulation scheme 
that subdivides the fiequaicy spectrum that it uses into a plurality of channels. The scheme 
^ves good performance over a dispersive transmission path. The data stream is split into 
Nso lower rate data streams that are transmitted over Nsd sub-carriers. The symbol duration 
increases for the parallel sub-carriers and hence the transmission scheme is more robust in 
the presence of multi-path interference. 

Finally the signal is transmitted by a radio frequency transmitter 16 such that the radio 
signal can propagate via a transmission channel 1 8 to a receiver 20. 



The channel may be dispersive and noisy. 



Consider the performance of a generic io, sub-carrier xised in the OFDM scheme. The ith 
channel carries a series of QAM symbols. If we consider just one of these, then 

a[i] = a,[q+jaQ[i] (1) 

Where 

a[i] represents a symbol in the ith sub-carrier 
a/[/] represents the real component of the symbol 
ag[/] represents the imaginary component of the symbol. 
j represents the imaginary operator. 

Each symbol represents a convolutionally encoded bit sequence given by 

6/.2 ; fti.m, bQXybgx. as shown in Figure 3 for a 16 QAM 

constellation.. 

The output of the receiver 20 is provided to an OFDM demodulator 22 whose output is fed 
to an equalizer and soft-output QAM demapper 24. Assuming that the cyclic prefix, 
introduced by the OFDM modulator 14, completely eliminates ISI (Inter OFDM Symbol 
Interference) and ICI (Inter Channel Interference) and that tiie channel estimate is error 
&ee» the received equalized signal in the id, sub-carrier is given by: 

Ai\ = r[i]/GcA(0 - a[i\ + H^[i]/Gc*(0 (2) 

Where 

r[i] is the received signal before equalization at the ou^ut of the OFDM demodulator 22 in 

Gch{i) is the Channel Frequency Response (CFR) coefficient (a complex number) in the ith 
sub-carrier, and 

is the complex Additive White Gaussian Noise (A WON). 
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The de-mapper 24 also receives information from a channel state information (CSI) 
estimator 26. The CSI estimator 26 attempts to deduce the effect of the transmission path 
on the channel. 

In broad terms the channel both scales the QAM signals and rotates them in phase ie the 
QAM constellation gets rotated. The CSI estimator attempts to estimate the effect of the 
channel on the symbol. The CSI estimator can woA on the assumption that channel is 
slowly time variant. However, the physical transmission scheme includes "preambles" 
which contain known sequences of symbols and these can be used to estimate the channel 
status. Thus individual phase and modulus estimates can be derived for each charmel on 
the basis of identifying the preamble. 

The output of the dem^er 24 is provided to a de-interleaver 28 and to a convolutional 
decoder 30. The Viterbi algorithm is a widely used method of decoding convolutional 
codes. The algorithm searches the possible code words of the convolutional code and 
detects the one that is most likely to have generated the received sequence. The search 
procedure steps through the code trellis and for each path along the trellis computes a 
metric which quantifies the discrepancy between the received sequence and the possible 
coded sequence. If the information associated to bits fed into the decoder is hard (ie 
binary, a sequence of -I's and Vs) the decoder is called hard decision decoder. 
Alternatively, if the information is soft, consisting of a hard decision (tiie sign) and a 
confidence level, or reliability (the magnitude), that represents how much confidence there 
should be in the hard-decision, the decoder is called soft-decision decoder. It's well known 
that, soft decision Viterbi decoding can give significant gain over hard decoding at the 
expense, however, of a greater computational complexity. In order to implement a 
soft-decision Viterbi decoder the demapper, which precedes the Viterbi decoder, needs to 
deliver soft information associated to the bits. 

Soft decision demapping for bit interleaved coded modulation (BICM) signals with BPSK 
or QPSK modulations is straightforward as the soft bit information, before being weighted 
by the Channel State Information (CSI) is simply given by the received signals for BPSK 
and by their in-phase and quadrature components for QPSK. Therefore, in the follovwng 



I 
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discussion consideration will be given to the higher modulation formats, for which soft 
detection requires much more computational effort. 



In the literature, two different approaches can be found to calculate the soft information for 
BICM signals, with multi-level modulations. 

The first prior art mechanism for bit interleaved coded modulation (BICM) schemes was 
disclosed by E. Zehavi, "8-PSK Trellis Codes for a Rayleigh Channel" IEEE Trans on 
Conun, Vol 40, pp 873 - 884, May 1992. The process starts by calculating sub-optimal bit 
metrics that are then used inside a Viterbi decoder for path metric computation. 

For each bit b\x and each bit bojc (where I and Q represent in-phase and quadrature parts, 
respectively, and K represents an index of the bit associated with the symbol where K is an 
integer in the range 1 ^ /l ^ /n) the QAM constellation is split into two partitions of 
complex symbols. 

These partitions are 

having symbols with 0 in position I, K. 
SjJ^ having symbols with 1 in the position I, K. 



Finally the metrics arc de-interleaved by a de-interleaver 28 and provided as an input to a 

Yiterbi d^oder 30. 

The Viterbi decoder works according to a well known algorithm which need not be 
described in detail here. However various web sites give tutorials in Viterbi decoding, 
such as http://Dwcb.netcom,com/-chip. f Adterbi/algrthms2.html. 




symbols with 0 in the position Q, K. 
symbols with I in the position Q, K. 




The bit metrics are given by 



A/c(6/,Jc)=ICc/K0l^ 



. min ly[ihal^ ♦ c=0,l 



(3) 



The convolutional encoder 2 functions as a state machine and the Viterbi decoder is 
furnished with a state map of the state machine showing which state to state transitions are 
allowed and which ones are disallowed. 

In the second prior art approach the QAM symbols are first demodulated by a soft ou^ut 
de-mapper and passed to a soft-input Viterbi decoder, see M. Speth et al, "Low Complexity 
Space-Frequency MLSE for Multi-User COFDM", IEEE GLOBECOM *99, pp 2395 - 99, 
Dec. 1999. 

In this approach the process seeks to de-map the received signal into soft bits which have 
the same sign as that provided by a hard decoder and whose magnitude indicates the 
reliability of the decision. 

The soft bit infonnation assigned to bit bijc can be shown to be given by the log-likelihood 
ratio (LLR) of the hard decision on bi.ic (see IL Pyndiah et al **Near Optimum Decoding of 
Product Codes", IEEE GLOBECOM 94, pp 339 - 43, Nov. - Dec, 1994) and can be 
approximated by 

LLRibf^K)^\GcHii)\^Drjc (5) 
We define 5' = 9i{sJl }as the subset containing the real parts of the complex symbols of 

(r) 



subset S^fl , for c = 0,1 . It can be shovm that equation (4) can be rewritten in a sunpler 



form. 



min (y/W-o/)^' min (y/W-fl/)^ I . 
ages'); '»/e5*); J 



Where the two minima are now taken over real values instead of complex ^mbols. 



If this is, for convenience, explicitly evaluated for the 16 QAM symbols we have 



>/W \yi[i]\<2 

^/j= 2(y/[z]-l yi[i]>2 

2(y;m + i) yi[i]<2 



Dia ="b/Wl+2 (7) 

equivalent expressions hold for the quadrature components with 'T' replaced by "Q". 

It has been demonstrated in F. Tosato and P. Bisaglia "Simplified Soft-Output Demapper 
for Binary Interleaved COFDM with AppUcation to fflPERLAN/2", IEEE ICC 2002, 
April-May, 2002. that using the approximate bit metrics Mc{bi^K) in equation 3 for path 
metric calculation inside the Viterbi algorithm is equivalent to demodulating the signals 
into soft bit values according to equation 6 and then employing a soft Viterbi algorithm for 
decoding. 

The formula for calculating the log likelihood ratio in equation 6 can be fiirther 
approximated by calculating |D/jc|(or indeed \Dqjc\) as the distance of the received 
equalised signal y{i\ from the nearest partition boundary within the partitioned QAM space 
and assigning Di^k (or Dq^k as appropriate) the sign or according to vrtiich partition 
y[i\ falls in. The magnitude (or absolute value) is a measure of distance of the received 
symbol from the decision boundary. 

Figure 4a illustrates the decision boundaries Bl to B4 for 16 QAM modulation and Figures 
4b to 4e illustrate the resulting partitions in 16 QAM space. The 16 QAM constellation is 
used by way of an example. However the present invention can be applied to higher order 
constellations in a similar way. 

Furthermore if we let dijc and dgjc denote half the distance between the partition 
boundaries B2 and B4 relative to bit bijc and bgx. respectively, then for the 16 and 64 

OAM cnn.^ellationK with fwrsiv nrtAnnino um\ in TPPPftno 1 1« ittdcdt axiai u« 

^ _ , - ^ J ^ - wa*> « « H A AAA A.««^<kM AA ak. W«AtA LFV 

shown that 

''«--|Z)„.,|**^.Jt>I 
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In terms of computational complexity the prior art system of decoding the symbols, even 
with ail of the simplifications and approximations invoked, is computationally complex. 

To illustrate this consider the case of a burst transmission. In such a burst transmission the 
channel can be assumed to be time-invariant for the duration of the burst. Thus channel 
state estimation need only be performed once by the receiver at the beginning of each burst. 

If we denote Nh the number of bits coded in the data burst and Nsd the ntmiber of 
sub-carriers {Ni, » Nsd) LLR calculation using the formulae: 

LLRib,jc) = \GcHii)\^Di,K 

LLRibQjd = \GcH(i)\^DQjc (9) 

requires one real multiplication per coded bit, plus computation of Nsd squared modulus 
of complex values per physical data burst, which is equivalent to real multiplications. 

Thus, the approximate LLR calculation requires {Nb + 2Nsd) real multiplications. 

According to a first aspect of the present invention there is provided a method of estimating 
a measure of trust of data conveyed by a QAM symbol^ wherein the measure of trust is 
calculated as a linear function of the modulus of a channel state estimation. 

It is thus possible to reduce the number of calculations, and in particular multiplications, 
performed in the soft output de-mapping by using an estimate of likelihood derived as a 
function of iGc/KOl rather than |Gc//(Ol^- 

Preferably the estimate of trust for a bit bjjc for K=I , where K represents a bit index within 
a complex symbol* is calculated as 1R{r[z]e^'''*^'^}, where 91 represents the "real'' part of 
a complex nimiber, r\l\ is the received symbol in the itb sub-carrier, and e'^^<*W represents 
the reciprocal of the phase response of an ith transmission channel over which the symbol 
was transmitted. 

Preferably the estimate of trust for a bit bgjc is calculated as 5{r[i]^"-'^'*^'^} where 3 
represents the imaginary part of a complex number, for K=l . 
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Preferably the estimate trust of a bit 6/^ is further calculated as 
-\LLRib,,K-i)\ +(IC7cft(0l djjc) for A" > 1 where Gc*(0 represents the channel frequency 
response complex coefficient on an iu, channel and d/^K denotes a half distance between 
partition boundaries in QAM space, for A' > 1, and the estimate of trust of a bit bgjc is 
calculated as -\LLRibQjc-i)\-^(\Gch[i]\*dQjc) where dg^ denotes the half distance 
between partition boundaries in QAM space, for A: > 1 . 

Preferably the estimate of trust is approximate log-likelihood ratio. 

According to a second aspect of the present invention, there is provided an apparatus for 
performing the method according to the first aspect of the present invention. 

The present invention will now further be described, by way of example only, with 
reference to the accompanying figures: 

Figure 1 schematically illustrates a transmit path and a receive path for a bit interleaved 
coded modulation scheme; 

Figure 2 schematically illustrates a convolutional encoder; 
Figure 3 illustrates a square 16 QAM constellation; 

Figure 4a illustrates the partition boundaries that are used to partition the 1 6 QAM space; 

Figures 4b to 4e illustrate the resulting partitions for a 16 QAM space; 

Figure 5 is a schematic diagram of a channel state estimator and soft output QAM 
demapper constituting an embodiment of the present invention; 

Figure 6 schematically illustrates a soft decision block constituting an embodiment of the 

nrftsent invpnHnn ff>r 1 ^ QAM COHSteU^tion; 

Figure 7 schematically illustrates a soft decision block constituting an embodiment of the 
present invention for 64 QAM; and 



( 
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Figure 8 is a comparison of decoding schemes, with the curves showing simulations for 
performance of a HIPERLAN/2 system. 



It should be noted that, for simplicity, the present invention is described with respect to 16 
QAM and 64 QAM, but it can be applied to larger constellations and similar demapping 
schemes can be derived. 

Unlike single carrier systems in which all symbols are affected by the same signal to noise 
ratio (on average), a multi-carrier OFDM system of the type shown in Figure 1 is such that 
each individual carrier suffers from an individual signal to noise ratio. However it is clear 
to the person skilled in the art that data conveyed on channels having a high signal to noise 
ratio is a priori more reliable than data transmitted on channels suffering from a low signal 
to noise ratio. This additional information has, in the prior art, been encoded by weighting 
the LLR functions by the square modulus of the channel frequency response, which 
represents the channel state information. 

However the inventors have realised that using the modulus of the channel fiequency 
response coefficients instead of the square of the modulus for calculating "soft bif * 
information for use by a decoder results in only a slight performance loss in terms of bit 
error rate at a given signal strengths as represented by Eb/No, where Eb is the energy per 
information bit and No is the power spectral density of the noise. 

This ^proximation does, however, allow for a computationally efficient implementation of 
a one tap equaliser and LLR calculation subsystem. Thus, the inventors have realised that 
sub-optimum soft input Viterbi decoding of a binary interleaved coded OFDM signal can 
be achieved with little additional complexity compared to realising the same operation 
using hard decoding instead. 

By applying this approximation to the prior art scheme for calculating the LLRs, the 
following equations are obtained. 

LLR(b .^^^'•W^"'''*^'^)' 

^ '-"^ -\LLR{b,Jc-^)\ + \GcH{i)\ -d/jc K>1 



n 



Where Gca(0 = ICc/,(Ole^^'*<'> 

Similar results hold for other Gray labelling patterns to that shown in Figure 3. 

Thus, compared with the prior art calculations of LLR, it can be seen that inside the LLR 
expression only the thresholds drjc and dgjc are scaled by the coefficients that convey the 
channel state information. 

As a consequence, the following scheme can be adopted for joint OFDM signal 
equalisation and LLR computation. A block diagram for implementing the scheme is 
illustrated in Figure 5. 

The incoming OFDM signal from a receiver is converted into individual data channels 1 to 
NsD by a Fourier transform block 50. The phase of the signal is then equalised/corrected. 
This is done by sending each one of the channels to an input of a channel frequency 
response estimator 52 which estimates the phase of each one of the channels and thereby 
produces a phase correction signal e for each channel i, where i is an integer in the 
range 1 ^ / ^ Nsd- The channel state estimate is done once per physical burst, at the 
beginning thereof 

The phase corrections are applied to each of the channels via respective multipliers 54(1) 
to 54(Nso). The phase equalised channel signals are then supplied to first inputs of 
respective soft bit calculators 56(1) to 56(Nsd). 

The channel frequency response estimator 52 also estimates the modulus of tiie signal 
iiieugui in each one of ihe channels, and this mlbrmation is supplied to a threshold 
calculator 58(1) to 58(Nsd) associated with each individual channel which calculates the 
threshold values lGc*(/)l • dijc and |Gca(i)I • dg,/: where i represents the channel number. 
These values are tiien passed to the respective soft bit calculators 56. 
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The specific implementation of the soft bit calculators is shown in Figure 6, for a 16 QAM 
constellation and in Figure 7 for a 64 QAM constellation. 

For the arrangement shown in Figure 6, the received phase equalised signal is received at 
an input 70 and supplied to a first analyser 72 v^^ich calculates the real component of the 
signal. The output of the first analyser represents bit value The value 

LLR(bij)is fiirther passed through an absolute value former 74 which calculates and 
negates the absolute value of LLR{bij). The ou^ut of the absolute value former 74 is 
then added to \Gch(/)\ • c// by an adder 76 to yield an output LLR(bi;i). A similar process 
is implemented on the imaginary component of the input signal to yield LL/2(6Q,i)and 
LLRibg;,). 

The arrangement shown in Figure 7 is similar to that shown in Figure 6 and like parts are 
denoted by like reference numerals. The additional features are that the output LLR{bi^ 
is made available to a fiirther adder absolute value former and negator 78 whose output is 
then added to \Gch{i)\ • ^2 by a further 80 to yield LLi?(fe 7,3). Similar components are 
provided in the "imaginary" signal path for calculating LLR{bi^)> 

Comparing the complexity of this scheme in a burst mode with the prior art, it is now 
necessary to calculate Nsd real multiplications per transmission burst in the case of 16 
QAM and 2Nso fo^ ^ QAM because only the thresholds dK need to be scaled by the 
cannel state information coefficients. 

However, these multiplications are not actual multiplications because for square QAM 
constellations dx are powers of two, so in fact only bit-wise shifts are needed for threshold 
scaling. 

It is, however, necessary to calculate Nsd moduli of the complex values per physical burst 
to calculate the channel state information values. This requires 2Nsd real midtiplications 
and Nsd square roots to be formed. 



Comparing the computational woridoad: 



13 

Present invention = 2Nso real multiplications -\-Nsd square roots. 
Prior art = INsd real multiplications -h^^ real multiplications. 

The computational advantage comes because Nh » Nso and in practice Nsd is 48 in 
HIPERLAN/2 whereas ^^is likely to be several thousand 

The performance of the present invention compared to the computationally complex prior 
art is shown in Figure 8. For a given bit error rate the degration in sensitivity is only a 
fraction of a dB for 16 QAM and is negligible for 64 QAM at BER 10"^. The 
approximation gets tighter for larger constellation size, where more calculations are saved 
by the approximate demapping. The gr^h compares hard decoding^ prior art soft decoding 
and soft decoding according to the present invention in modes 5 and 7 of HIPERLAN/2 
using 16 and 64 QAM respectively. 
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CLAIMS 

A method of estimating a measure of trust of data conveyed by a QAM symbol, 
wherein the measure of trust is calculated as a linear ftmction of the modulus a 
channel state estimation. 

A method as claimed in claim 1, wherein the measure of trust is a log-likelihood 
ratio. 

A method as claimed in claim 2 where the ILJ?(&y^), for K=l, is calculated 
as 91 {r[/]e~^^^^}, where 91 represents the "real" part of a complex number, r[i] is 
the received symbol in the id, channel, and e'^^^^ represents the reciprocal of the 
phase response of the itb transmission chaimel over which the symbol was 
transmitted, where K represents a bit index within a complex symbol. 

A method as claimed in claim 3, where the LLR(bQjc) for K=l is calculated as 
3{r[z]e''^**^'^} where 3 represents the imaginary part of a complex number. 

A method as claimed in claim 3, where the LLR{bix) is calculated as 
-|LLiJ(6/jc-0l +(lGcir(Ol •^'/jc) for *^>1 where Gca(0 represents the channel 
frequency response complex coefficient on an ia channel and di/c denotes a half 
distance between partition boundaries in QAM constellation. 

A method as claimed in claim 3, where the LLRQygjc) is calculated as 
'-\LLR{bQjc-\)\^(i^Gch{i\\'dQjc) for K>\ where dgjc denotes the half distance 
between partition boimdaries in QAM constellation. 

A method as claimed in claim 3, wherein thresholds dfjc are scaled by the amplitude 
of the channel state estimator coefficients Gca(0* 

A method as claimed in claim 7, in which computational complexity is reduced by 
using the modulus of the estimated channel state estimator coefficients IGcaO)! 
instead of |Gca(OI^ ^ channel state information inside a soft output demapper. 

A method as claimed in claim I, wherein the QAM symbol is a symbol in a OFDM 
transmission scheme using multi-level modulation. 
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soft output demapper operating in accordance with the method claimed in claim 
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